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Abstract: We obtain several new results for the complex generalized associated Lame potential 
V{x) = a{a + l)m sn^(y, m) + b{b + l)m sn^(y + K{m),m) 

+f{f + sv?{y + K(m) + iK'{ni),m) + g{g + l)m sn^(?/ + iK' {m),m) , 

where y = x — ^^^^ — i^JlI!) ^ sn(y, m) is a Jacobi elliptic function with modulus parameter m, and there 
are four real parameters a, b, f, g. First, we derive two new duality relations which, when coupled with 
a previously obtained duality relation, permit us to relate the band edge eigenstates of the 24 potentials 
obtained by permutations of the four parameters o, b, f, g. Second, we pose and answer the question: how 
many independent potentials are there with a finite number "a" of band gaps when o, 6, /, g are integers? 
For these potentials, we clarify the nature of the band edge eigenfunctions. We also obtain several analytic 
results when at least one of the four parameters is a half-integer. As a by-product, we also obtain new 
solutions of Heun's differential equation. 



1 Introduction. 

In a recent paper PP, hereafter referred to as I, we discussed the generalized associated Lame (GAL) 
potential given by 

V{x) = a{a + l)m sn^ {x , m) + b{b + l)m sn^ {x + K (m) , m) 

+ /(/ + l)m sn^(a; + K{m) + iK' {m),m) + g{g + l)m sn^(a; + iK'(m),m) 

/ -,^ 9/ N ,/, -,^ cn^(2;,m) ^ , ^ ,dn^(x,m) , 1 

= a{a + l)m sn2(x, m) + 6(6 + l)m ^) [ + f{f + 1) ( + g{g + I)—, r , 

dn (x,m) cn^(a;,mj sn^(x,mj 

(1) 

which involves four real parameters a,b,f,g. Here, sn {x,m), cn (x,m), dn (x,m) are Jacobi elliptic 
functions with elliptic modulus parameter m (0 < m < 1). They are doubly periodic functions with 
periods [4:K{m),i2K'{m)], [iK (m) , 2K (m) +i2K'{m)], [2K{m),i4K'{m)] respectively P], where K{m) = 
Jq^'^ d6[l — msin^^]^^/^ denotes the complete elliptic integral of the first kind, and K'{m) = K{1 — m). 
Prom now on, unless essential, we will not explicitly display the modulus parameter m as an argument of 
Jacobi elliptic functions. It may be noted here that the four terms in the GAL potential correspond 
to complex translations of the independent variable x by 0, K{m), K(m) + iK'{m),iK'{m). Although the 
GAL potential is real, it does have singularities on the real axis coming from the zeros of the Jacobi elliptic 
functions sn (x) and cn (x) in the last two terms. One way to avoid these singularities is to make a complex 
change of variables y = ix -\- (3, with [5 being an arbitrary real, non-zero constant, chosen so as to avoid the 
singularities arising from the zeros of Jacobi elliptic functions on the real axis ,3^ . This procedure was used 
in I, and we studied many of the properties of the resulting PT-invariant complex periodic potential 1^. 
However, consistent with the practice in the mathematics literature, in this paper we use an alternative 
approach of avoiding singularities by simply translating the independent variable x by an arbitrary non-zero 
amount in the complex plane. In fact, for simplicity, from symmetry considerations, we take the translated 
variable to be x — ^'^^^ — , Note that the potential still is a PT-invariant complex periodic potential. 

One very important consequence of this choice is that the energy eigenvalues which we shall obtain here 
will be opposite in sign from the values obtained in I. This point should be kept in mind while comparing 
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the results from the two papers. Thus exphcitly, we consider the potential 



V{x) 



y{y) 



a{a + l)m sn^(y, m) + b{b + l)m sn^(y + K{m),m) 




sn- 



1 



[a, b, f, g] 



(2) 



where 



y = x- 



K{m) iK'{m) 
~ 2 



(3) 



It may be noted here that in I we used the notation [a(a + 1), b{b + 1), /(/ + l),g{g + 1)] to denote this 
potential. However, for consistency with the prevailing practice in the mathematics literature, in this paper 
we use the notation [a,b,f,g]. In this notation, ordinary Lame potentials are denoted by [a, 0,0,0], and 
associated Lame (AL) potentials are denoted by [a,6, 0, 0]. Note that the potential ^ remains unchanged 
when any one or more of the parameters a, b, /, g change to —a — 1,—b — 1, — / — 1, —g — 1 respectively. 

There is one important point which permits us to construct many supersymmetric partner potentials 
corresponding to a given potential [31 . Although this point was previously made in I, it is worth restating 
here. Normally, in supersymmetric quantum mechanics 6 , given a potential V-{x), the ground state 
wave function V'o(a^) is used to construct the superpotential W{x) = —■iP'q(x)/iPo{x), which then yields the 
supersymmetric (SUSY) partner potential V"+(x) = W'^ + W. If one uses any excited state wave function 
■0(2;) of V-{x) to construct a superpotential W{x), then the original potential V^{x) is recovered correctly 
(by construction), but the corresponding partner potential V+{x) turns out to be singular on the real x- 
axis due to the zeros of the excited state wave function tp{x). However, for the complex potential ((2), the 
singularities are not on the real axis L7J. Thus in this case one could also use any of the excited state wave 
functions to obtain the superpotential W{x) and hence discover several supersymmetric partner potentials 
with the same energy spectrum. 

In I we showed that several GAL potentials with specific integer values of a, b, f, g have a finite number 
of band gaps. Further, looking at the symmetry of these potentials, we conjectured that all GAL potentials 
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with integer values of a, b, f, g also have a finite number of band gaps. Some results of this type are available 
in the mathematics literature, and it is worthwhile to present a brief review of what is already known about 
the GAL potential. 

The GAL potential @, expressed in terms of Weierstrass functions, was discussed in a brief note by 
Darboux in 1882 8 , as well as in two subsequent articles in 1914 and 1915. He mentioned that some 
results had already been presented by Hermite in 1872 in unpublished lectures at Ecole Polytechnique. In 
1883, Sparre wrote two long papers on the GAL potential expressed in terms of Jacobi elliptic functions. 
Unfortunately, the mathematics community was largely unaware of these papers, and the GAL potential 
was "rediscovered" by Treibich and Verdier in 1990 jJUj. In the current mathematics literature, the GAL 
potential, expressed in terms of Weierstrass functions, is known as the Treibich- Verdier potential. Several 
workers have shown that when a,b,f,g are integers, these potentials have a finite number of band gaps. 
In particular, the number of band gaps p is given by 



1 

p = — max 
^ 2 



2 max [a, 6, /, ff], 1 + TV - 1 + (-1)^^ min [a, b, f,g] + - 



1 



(4) 



where N = a + b + f + g. It was also shown that finite band gap potentials are solutions of higher order 
KdV equations 12 . In recent times, several people have also discussed the connection between Heun's 
equation and the Treibich- Verdier potential |l,Sj . Finally, some authors have studied even more general 
potentials with a finite number of band gaps jl4j . 

In I using supersymmetry we showed that the band edge eigenvalues of the Lame and the AL potentials 
[2, 0, 0, 0], [3, 0, 0, 0], [a, (a — 3), 0, 0] are the same as those of the GAL potentials [1, 1, 1,0], [2, 1, 1, 1], [a, a — 
1,1,0] respectively. We had also conjectured in I that for integer a,b, Lame and AL potentials have the 
same band edge eigenvalues as some GAL potentials - the explicit relationship being: 

'a a a a — 2' 



[a, 0,0,0] = 
[a, 0,0,0] = 
[a, 6, 0,0] = 
[a, 6, 0,0] = 



2'2'2' 2 
a + 1 a — 1 a — 1 a — 1 



a = even integer 



2 ' 2 ' 2 ' 2 
a + b a + b a — b a — b — 2 



2 ' 2 ' 2 ' 2 
a + b + 1 a + b—1 a — b 



1 a 



a = odd integer , 

, a + b = even integer , 

a + b = odd integer . (5) 



1 
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Note that we are using the notation "=" to denote "same band edge eigenvalues", but not identical 
potentials. It is interesting to note that recently, Takemura ^15: has verified the conjectures expressed in 
©. Further, more generally, he has proved that the GAL potential [a,b,f,g] for integer a,b,f,g has the 
same band edge eigenvalues as another GAL potential, with the explicit relationship depending on whether 
N = a + b + f + g is an even or an odd integer. If N is an even integer, the relationship is 

'a + b + f-g a + b + g-f a + f + g-b b + f + g-a 



[a, b, /, g] 

while if is an odd integer, then the relationship is 

'a + b + f + g + 1 a + b — f — g — 1 a + f — b — g — 1 a + g 



(6) 



[«, f, g] 



b-f-i 



(7) 



He has also shown that if a, b, f, g are all half-integers and their sum is an even integer, then the band edge 
energy eigenvalues of the potential [a, b, f, g] are the same as two other GAL potentials where also all four 
parameters are half-integers with their sum being an even integer. In particular, the explicit relationships 
are 



k -\ — ,6 = /H — , f = n -\ — ,g = »H — 
2 2 2 2 



k + l + n + p + ?, k + l 



n 



p—1 k + n — l—p—1 k+p- 



n 



l-l 



+l+n-p+l k+l 



2 ' 2 ' 2 

■p — n + 1 k + n + p—l + 1 I -\- n + p — k -\- 1 



(8) 

The above mentioned results, when combined with our work in I, raise several questions. For example, 
while it is clear from the relation @ that all the 24 potentials obtained by permutation of the four 
parameters a, 6, /, g have the same number of band gaps, what is the precise relationship between the band 
edge eigenstates of these 24 potentials? Second, a very interesting question is to ask how many independent 
GAL potentials are there with say a band gaps. Third, is it possible to further generalize the results of 
Takemura |15j ? Besides, what is the nature of the band edge eigenfunctions for a general GAL potential 
in case a, 5, /, g are all integers? Finally, in view of the connection between the GAL potentials and Heun's 
equation [^1 , it is worth enquiring about the implications of these results in the context of solutions to 
Heun's equation. 
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In this paper, we address all the above-raised issues. Further, we also consider GAL potentials in which 
at least one of the four parameters is a half-integer while the other parameters are arbitrary numbers 
and show that both relations © and O are valid in all these cases. As a consequence, we conjecture 
that the relations © and Q are in fact simultaneously valid even when all four parameters a, b, /, g are 
integers with their band edge eigenvalues being the same as the quasi-exactly solvable (QES) eigenvalues 
of potentials where all four parameters are half-integers (with their sum being an odd integer). Further, 
knowing the QES eigenvalues of the GAL potentials [a, 1/2, /, g], [a, 6, 1/2, g], [a, 6, /, 1/2] (where a,b,f,g 
are arbitrary numbers), and using the connection between the Schrodinger equation for the GAL potentials 
and Heun's equation, we show that the corresponding eigenfunctions can be obtained by directly solving 
the algebraic form of Heun's equation. 

The plan of this paper is as follows. In Sec. 2, we derive new duality relations and examine some 
consequences. In particular, using these duality relations we obtain the precise connection between the 
energy eigenvalues and eigenfunctions of the 24 potentials obtained by permuting the four parameters 
a,b,f,g. In Sec. 3, we discuss the GAL potential @ in some detail when the four parameters a,b,f,g 
are all integers. First, we find various GAL potentials which are related to each other in the sense that 
they have identical band edge energy eigenvalues. We also find a large number of self-isospectral potentials 
as well as self-dual but non-self-isospectral potentials. Using all these results, we obtain the number of 
independent potentials with say a band gaps. We also clarify the nature of the band edge eigenfunctions 
for these potentials. In Sec. 4 we discuss the GAL potential when either one or more of the four parameters 
take half-integer values while the remaining parameters are arbitrary. In this case, in general one expects 
to obtain QES mid-band states, by which we mean any energy state lying inside an energy band and thus 
not a band edge. Quite remarkably, by generalizing Takemura's results |15| we obtain GAL potentials 
which have the same band edge eigenvalues as the mid-band energy values of the potentials with one or 
more of the parameters being half- integers. In Sec. 5 we discuss the GAL potential in case the parameters 
a, b, /, g take arbitrary values and obtain the corresponding GAL potentials with the same band edge and 
mid-band energy values when a + b + f + g [or any other combination obtained by replacing one or more 
of these parameters by —a — 1,—b — 1, — / — 1, —g — 1 respectively] is an even integer. In Sec. 6 we discuss 
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the implications of all these results in the context of Heun's equation. In particular, we show that in 
view of the connection between the different GAL potentials, given a periodic solution of Heun's equation, 
one immediately obtains four periodic and three quasi-periodic solutions of the same equation. We also 
show that in many cases, knowing the QES energy values of a GAL potential, it is much easier to solve 
the algebraic form of Heun's equation and obtain the corresponding eigenfunctions of the GAL potential. 
Finally, in Sec. 7 we summarize the results obtained in this paper and spell out some open problems. 

2 Duality Relations for GAL Potentials. 

In this paper our main focus is on the Schrodinger equation 

- —.i^ix) + V{x)'iIj{x) = EiP{x) , (9) 

where V{x) is the GAL potential given by eq. and we have chosen units with h = 2m = 1. Displaying 
parameters more explicitly, eq. © states that the potential V{a, b, /, g, m; x) has band edge eigenvalues 
E{a,b, f, g,m) and eigenfunctions 'ip{a,b, f,g,m;x). Of course, this means that the translated potential 
V{x) = V{a,b, f,g,m;x) = V{a,b, f, g,m;x — — i^JH-j given in eq. Q has the same eigenvalues 
E{a, b, f, g, m) and translated eigenfunctions ip{a, b, f, g,m;x — — lIL^Ill'j^ 

First we want to show that the band edge eigenstates of the 24 potentials [obtained via permutations 
of the 4 parameters, a,b,f,g in V{x)] are all related, so that once the band edge eigenstates of any one 
permutation are known, the complete band edge eigenstates of all 24 potentials are also known. Actually, 
these relations are valid for both band edges as well as mid-band states. 

The first relation is simple - in view of the invariance of the Schrodinger eq. © under the translation 
x ^ X + K{m) followed by the interchanges a ^ b and f g, one gets P 

E{b,a,g,f,m)=E{a,bJ,g,m), ip{b, a, g, f,m; x) oc i;{a,bJ,g,m;x + K{m)). (10) 

Similarly, the translations x ^ x + iK'[m) and x — > x -|- K{m) + iK'[m) followed by suitable interchanges 
of parameters, yield 

E{g,f,b,a,m)=E{a,b,f,g,m), Tp{g, f,b,a,m;x) cc ip{a,b, f, g,m;x + iK'{m)) , (11) 
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E{f,g,a,b,m)=E{a,b,f,g,m), tp{f,g,a,b,m;x) oc ip{a,b, f, g,m; x + K{m) + iK'{m)) . (12) 

Thus, once we obtain the eigenvalues and eigenfunctions of a given potential [o, b, f, g], then we immedi- 
ately know the eigenvalues and eigenfunctions of three other potentials: [b, a, g, /], [g, f, b, a] and [/, g, a, b]. 
Note that relations (|1UI) . all involve the same modulus parameter m. 

We now derive three remarkable duality relations, which connect the energy states of different GAL 
potentials, and involve changes in the modulus parameter from m to 1 — m, 1/ni and — — m). 
Duality Relation I: This was already derived in I and is given by 

Eia, b, f, g, m) = [a{a + 1) + 6(6 + 1) + /(/ + 1) + g{g + 1)] - E{a, g, /, 6, 1 - m) , 

ip{a,b, f, g,m;x) oc ip{a, g, f,b,l — m;ix + K'{m) + iK{m)) . (13) 

Duality Relation II: Using the formulas |2] 

sn(x,m) = —;=sn{^/rnx, — ) , cn(2;,m) = dn{^/m.x, — ) , dn(x,m) = cn{^/mx, — ) , (14) 



mm m m 



and redefining a new variable z = -y/mx, the Schrodinger eq. Q takes the form 



m mm dn^(2;, ;i^) cn^(z, ^) 



, , , E(m) , , , , 

Hz) = -^i^iz) . 15 
m 



sn^(z, — ) 

On comparing eqs. © and (fT3|) we obtain the duality relation 

E{a,b,f,g,m) = mE{a,f,b,g, — ) , xlj{a,b, f, g,m;x) oc ^jia, f,b, g, —■,^/mx) . (16) 

m m 

Duality Relation III: We again start from the Schrodinger equation Q and now use the formulas |2j 



sn[-v/l — m x, jz 

sn(x, m) — 



\/l — m dn[vT" 



m X 



' l~m' 



cn[Vl - m X, j^] 

cn(x,m) — 



dn[\/l — m X, 



l—m' 

dn(x,m) = 3—. (17) 

dn[Vl - m X, jr^l 
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On defining a new variable z = \/l — m x, the Schrodinger eq. © takes the form 

+ /(/ + ! 



(f ^, . . \ m 6(6 + 1) 2/ -™ N m o(a + l) cn2(z,^) ^ ^ ^^dn2(z,^) 



1 — m 



-sn (z, 



1 — m 1 — m dn (z, j^^) 



+ 9(9 + 1) 



' 1— m ' 



1 



a(a + 1) + 6(6 + 1) + /(/ + 1) + gig + 1) 



1 — m ^ ^ ' 1 — m 
Comparing eqs. © and (fT8|) . one gets the duality relation 



^(z). 



(18) 



£;(a, 6, /, 5, m) = (1 - m)£;(6, a, f, g, 



-m 



1 — m' 



+ m 



a{a + 1) + 6(6 + 1) + /(/ + 1) + g{g + 1) 



^'(a, /, 5> ni; x) oc ?/;(6, a, /, g, 



Vl — "I x) 



(19) 



1 — m 

Using the three duality relations [eqs. (|T^ . ifTB]) . (fT^ ] along with the translation results [eqs. (fTU]) . ((TT|) . 

(|12p]. it is easily shown that once the eigenstates of a given potential are known, we can immediately obtain 
the energy eigenstates of all the 24 potentials obtained by permuting the 4 parameters a, 6, /, g. Hence, out 
of these 24 potentials, there is only one independent potential and without loss of any generality, throughout 
this paper we only consider the potential [a, 6, f,g\ with a> b > f > g (unless stated otherwise). 

The duality relations are very powerful and have many interesting consequences. For example, we find 
that for arbitrary integer values of a, f, the potential [a, 0, /, 0] has only a finite number of band gaps. This 
follows because, eq. (|TB]) gives 



E{a,0,f,0,m) =mE{a,f, 0,0,-) 

m 



(20) 



so that both the potentials [a, /, 0, 0] and [a, 0, /, 0] must have the same number of band edges. Since one 
knows |16j that the AL potential [a, f, 0, 0] has a finite number of band gaps, the same statement holds for 
the potential [a, 0, /, 0]. 



3 Independent GAL Potentials with a Band Gaps [a,b,f,g = integers]. 

In this section, we want to answer the following interesting question: given a potential of the form [a, 6, /, g] 
with a>b>f>g>0, how many independent GAL potentials are there with exactly a band gaps? 
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To begin, for a given integer a, let us calculate the total number of possible GAL potentials with 
a>b>f>g>0. The number of such potentials is 

EEEl = (« + l)(« + 2)(« + 3)/6. (21) 

6=0 /=0 3=0 

Now due to the Landen transformations 121 a potential [a, a, b, b] with a 7^ 6 is essentially the same 
as the potential [a, b, 0, 0] and hence not really distinct. There are a such potentials. Similarly, the potential 
[a, a, a, a] is related to the potential [a, 0, 0, 0] and hence not distinct. This follows from the relations |2] 

sn {x) = ; — - — - , msn (x + %K {m),m) = — ^ , (22) 

^ ^ l + dn(22;) ' ^ v 7 sn2(x,m) ' ^ ^ 

sn'^{x,m)+sn'^{x+K{m),m)+sn'^{x+iK\m),m)+sn'^{x+K{m)+iK'{'m),m) = Asn'^{2x+iK'{m)) . (23) 
Thus the number of distinct GAL potentials of the form [a, b, /, g] with a>6>/>(7>Ois given by 

Ndist = a{a + l){a + 5)/6 . (24) 

In view of the formula Q for the number of band gaps, it is obvious that these distinct potentials are 
of two types - there are those with exactly a band gaps and those with more than a band gaps. Using eqs. 
© and Q) it follows that in case a + b + f + g is an even integer and a + g > b + f, then the potential has 
a band gaps, while if a + g < b + f, then it has {a + b + f — g) /2 band gaps. Similarly, using eqs. Q and 
(llj it follows that in case a + 6 + / + is an odd integer and a > b + f + g, then the potential has a band 
gaps, while if a < b + g + f, then it has {a + b + f + g + l)/2 band gaps. 

It is now straightforward to count the number of independent GAL potentials [a, b, f, g] with a band 
gaps and show that 

Na = ^[a^ + 9a^ + Qa + 2] , a = l (mod 3) ; 
18 

N^ = —[a^ + ga^ + Qa-2] , a = 2 (mod 3) ; 
18 

iV^ = -Ji[a2 + 9a + 6] , a = (mod 3) . (25) 
18 

Several comments are in order at this stage. 

1. There is only 1 independent potential with one band gap but there are 3 independent potentials with 
2 band gaps and 7 independent potentials with 3 band gaps. This implies ^2] that while there is only 
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one independent KdV equation of 3rd order, there should be 3 such equations of 5th order and 7 such 
equations of 7th order. It is worth pointing out that indeed there are 3 independent equations of 5th 
order cahed Sawada-Kotera and Kaup-Kupershmidt equations. It may be interesting to exphcitly 
obtain all the 7 independent KdV equations of 7th order. 

2. The number of independent potentials of the form [a, b, f, g] having more than a band gaps is obtained 
by subtracting eq. (^5]) ) from eq. H24() : 

7V>a = — (a + 2)(2a2 + 5a-l), a = 1 (mod 3) ; 
18 

N>a = ^{a + l)(2a2 + 7a + 2), a = 2 (mod 3) ; 
lo 

i\r>a = -3-a(a + 3)(2a + 3), a = (mod 3) . (26) 
18 

For example, for a = 2, there are 4 potentials of the form [a, b, f, g] with more than 2 band gaps: 

Ar>2 = [2,2,2,0], [2,2,2,1], [2,2,1,0], [2,1,1,1] . (27) 
For a = 3, there are 9 potentials with more than 3 band gaps: 

iV>3 = [3,3,3,0], [3,3,3,1], [3,3,3,2], [3,3,2,0], [3,3,2,1], [3,3,1,0], [3, 2, 1, 1], [3, 2, 2, 0], [3, 2, 2, 2] . 

(28) 

3. If a + 6 + / + 5 is an even integer and further if a + g = 6 + /, then both sides of eq. © are 
identical and the corresponding potential is self-dual. Similarly, if a + 6 + / + 17 is an odd integer 
and a = b + f + g + 1, then both sides of eq. are identical, and one has a self-dual potential. 
As an illustration, the AL potential [a, a — 1,0,0] is self-dual, a fact we already knew from ref. jl8j . 
But we now get a large number of additional self-dual GAL potentials, like [2, 1, 1, 0], [4, 2, 2, 0], and 
[4, 2, 1, 0], for example. 

4. A related question is, out of the above self-dual potentials how many are also self-isospectral [TSIIT^ ? 
It is easy to see that a choice of eigenfunction yielding self-isospectral potentials is 

dn"(x, m)sn^(a;, m) 
^ = b/ \ ' (29) 

11 



and that the corresponding self-isospectral potential has the form 

[a,a-l,b,b-l] , b > ; [a, a -1,0,0], 6 = 0. (30) 

The fact that the AL potentials of the form [a, a — 1,0,0] are self-isospectral was established many 
years ago [21 • However, what is new is the realization that the potentials [a, a — l,b,b — 1] are also 
self-isospectral for arbitrary integer values of a, b. Some examples are [1, 0, 0, 0] , [2, 1, 0, 0] , [3, 2, 2, 0] , 
[3,2,1,1], and [3,2,0,0]. Thus out of all the self-dual potentials, those which are of form (|30() are 
also self-isospectral, while the rest are self-dual but not self-isospectral. The number of self-dual 
potentials A^^^^ which are not self-isospectral is 

Nsd = ^{a-lf , a = 1,4, 7,... ; 

Nsd = ^a{a-2) , a = 2, 3, 5, 6, ... , (31) 

whereas the number of self-isospectral potentials is Nsi = a. As an illustration, the self-isospectral 
potentials with two band gaps are [2,1,0,0] and [2,1,1,0] while there are no self-dual, non-self- 
isospectral potentials with two band gaps. On the other hand, the self-isospectral potentials with 
three band gaps are [3,2,2,1], [3,2,1,0], and [3,2,0,0], while the only self-dual, non-self-isospectral 
potential with three band gaps is [3, 1, 1, 0]. 

5. Clearly all potentials with a + g>b + fora>b + f + g + l depending on if a -|- 6 -|- / -|- 51 is an odd 
or an even integer, have partner potentials as given by eqs. © or (jT)) respectively. Now out of these, 
some are SUSY partner potentials while the rest are merely partner potentials. So let us count both 
types of potentials. Now if two GAL potentials are SUSY partners, then one of their eigenfunctions 
must be related to each other by ijjn = i/jj^ ■ Further, these two eigenfunctions must have the form 

dn°(x)cn^(x)sn^(x) . (32) 

For such eigenfunctions, we know from I that b = —a,f = —(3,g = —7 and a + b + f + g = 0. 
From here it is easy to show that the potential [a, 6, /, g] with a band gaps has a SUSY GAL partner 
provided either a + g = b + f + 2 or a = b + f + g + Z depending on if a -|- 6 + / -|- is an even 
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or an odd integer respectively. Hence, all potentials of the form [a, b, f, g] with a band gaps and 
satisfying a + g > 5 + / + 2ora > b + g + f + 3 have merely partner potentials of the form © 
or (0) respectively. It is worth emphasizing that while these partner potentials have the same band 
edge eigenvalues, none of them are SUSY partner potentials. For example, [2,0,0,0] and [1,1,1,0] 
are SUSY partner GAL potentials with two band gaps. Similarly [3, 1, 0, 0] and [2, 2, 1, 0] are SUSY 
partner GAL potentials with three band gaps. On the other hand, [4, 0, 0, 0] and [2, 2, 2, 1] are merely 
partner potentials with four band gaps. 

6. One can count the number of potentials {Ngu) of the form [a, b, f, g] with a band gaps having another 
GAL potential as its SUSY partner and it is easy to show that 

N,^ = ^[a^-1] , a/0 (mods) ; 

Nsu = \a^ , a = (mod 3) . (33) 

Finally, it is not difficult to show that the number of potentials of the form [a,b,f,g] with a band 
gaps and having merely a (non-SUSY) partner potential of the form as given by eqs. © and ((T)) 
respectively is given by 

Nnsu = ^{a - l){a^ - 2a - 2) , a = 1 (mod 3) ; 

io 

Nnsu = 7^(« - 2)(a^ - a - 2) , a = 2 (mod 3) ; 

io 

Nnsu = ^a^{a-3) , a = (mod 3) . (34) 

io 

7. On adding the number of self-dual (but non-self-isospectral), self-isospectral, SUSY partner and (non- 
SUSY) merely partner potentials as given by eqs. (|5n|) to respectively, as expected, we find that 
the number of independent potentials with a band gaps is as given by eq. H25() . 

8. One obvious interesting question is whether there are non-GAL potentials with a finite number 
of band gaps. The answer to this question is yes. In particular, since the general form of the 
eigenfunction for any GAL potential is of the form 

N 

TpGALix) = dir''{x)cn~^{x)sn~'^{x) ^ Aksn^''{x) , (35) 

k=0 
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it follows that non-GAL potentials of the form 

V+{x) = Vgal{x) - 2—r \x^^gal{x) , (36) 

are also finite gap potentials. In this context it is worth mentioning that some people '14] have 
recently obtained potentials with a finite number of band gaps which are more general than the GAL 
potentials. It is not clear if those potentials and the potentials (|36() have any overlap. 

3.1 The Nature of Band Edge Eigenstates 

In I we showed that iia + b + f + g = 2n, then n + 1 QES states can be obtained for GAL potentials 
and they are of the form given in eq. (|35j) . On using this key result as well as the fact that the GAL 
potential [a, b, f, g] remains unchanged when any one (or more) of the four parameters a, b, f, g is changed 
to —a — l,—b — l,—f — l,—g — l respectively, it is easy to specify the nature of the band edge eigenfunctions 
for any potential [a, b, f, g] with a band gaps. The nature as well as the number of eigenstates crucially 
depend on whether a + b + f + g is an even or an odd integer - so we will discuss these situations separately. 
a + b + f + g = even integer: 

In this case it is easy to show that one has {a + b + f + g + 2) /2 eigenstates of the form 

sn"^^(2;)cn~^(x)dn"''(x)F(a+fc+j+g)/2[sn2(x)] , (37) 
(a + b — f — g) /2 eigenstates of the form 

sn^+i(x)cn^+i(x)dn-^(x)F(„+6-/-3-2)/2[sn2(x)] , (38) 
(a + f — b — g)/2 eigenstates of the form 

snS+\x)cn-f{x)dn'+Hx)F^,+f_,_g_,y,[sn\x)] , (39) 
and (a + g — b — f)/2 eigenstates of the form 

sn-5(x)cn^+i(x)dn^+i(x)F(,+g_,_;_2)/2[sn2(x)] . (40) 
If instead b + g > a + f , then one has {b + f — a — g)/2 eigenstates of the form 

sn^+\x)cn-f{x)dn~'{x)F^,+f_a-g-2)/2[sn\x)] . (41) 
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a + b + f + g = odd integer: 

In this case it is easy to show that one has {a + b + f — g + l)/2 eigenstates of the form 

sn3+\x)cn-f{x)dn-\x)F^^+,+f_g_^y2[sn\x)] , (42) 
(a + b + g — f + 1) /2 eigenstates of the form 

sn-9{x)cnf+\x)dn-\x)F^a+b+g-f-i)/2[sn\x)] , (43) 
(a + f + g — b+l)/2 eigenstates of the form 

sn-^'(x)cn-/(x)dn^+i(x)F(,+^+,_,_i)/2[sn2(x)] , (44) 
{a — b — f — g — 1) /2 eigenstates of the form 

sn^+i(x)cn/+i(x)dn^+i(x)F(,_,_^_,_3)/2[sn2(x)] , (45) 
If instead b + f + g>a— 1 then one has (b + f + g — a + 1) /2 eigenstates of the form 

sn-3{x)cn-f{x)dn-\x)F^b+f+g^a^i)/2[sn\x)] , (46) 

Here -F„[sn^(x)] denotes a polynomial of order n in sn^(x). It is worth pointing out that not only the QES 
band edge eigenvalues are identical for the two partner potentials as given either by eq. ® or by 0, even 
the nature of the band edge eigenfunctions in the two cases is also similar. For example, the potential 
[{a + b + f-g)/2, {a + b + g- f)/2, {a + f +g-b)/2, {b + f + g-a)/2] as given by eq. © has the same band 
edge eigenvalues and further the corresponding eigenfunction is simply obtained from eqs. (|37() to ()41() by 
replacing a, b, f, g with {a + b + f — g)/2,{a + b + g — f )/2,{a + f + g — b)/2,{b + f + g — a)/2 respectively. 
Exactly the same is also true about the equivalent potential given by eq. Q in case a + b + f + g \s w. 
odd integer and the corresponding eigenfunctions are exactly as given by eqs. (|42() to (|46() but with the 
replacement of a, 6, /,g by {a + b + f + g + l) /2,{a + b- f - g -I) /2,{a + f -b- g -I) /2,{a + g-b- f -I) /2 
respectively. It follows from here that irrespective of whether a + b + f + g \s w. odd or an even integer, 
there are precisely a bound bands, same a number of band gaps and 2a + 1 number of band edges all 
of which are analytically known in principle, beyond which there is a continuum band extending up to 
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E = oo. Further, irrespective of whether a + b + f + g is an odd or an even integer, if a + 6 is an even (odd) 
integer, then there are a + b + 1 band edges of period 2K{4:K) and a — h band edges of period AK{2K). 

Thus in general the band structure of the GAL potentials is unusual in that if a + 6 is an even (odd) 
integer, then h band gaps of period 4:K{2K) must be of zero width, i.e. there must be h doubly degenerate 
states of period 4K{2K). Unfortunately, till today we do not know either the eigenvalue or the nature of 
the eigenfunction of even one of these doubly degenerate states. One exception is the case of pure Lame 
(and their GAL partners) potentials, i.e. when b = f = g = 0, as in that case depending on if a is even 
or odd integer, one has a + 1 band edges of period 2K{4K) and a band edges of period AK{2K) and the 
band structure is normal one, with no doubly degenerate states. 

As an illustration, consider the GAL potentials with two band gaps. As seen above, there are 3 distinct 
potentials with 2 band gaps out of which we have already discussed the band structure of the two potentials 
[2, 0, 0, 0] and [2, 1, 0, 0] |16lll8j . Thus it would be interesting to know the band edges and the band structure 
of the remaining potential with two band gaps, i.e. [2, 1, 1,0]. Since a + 6 = 3, it follows from the above 
discussion that in this case there must be 4 band edges of period AK and 1 band edge of period 2K. Using 
Table 4 of I it is easily seen that the eigenstate with period 2K is given by 

V' = dn^(2;)sn(x)cn~^(x) , E = 9m , (47) 

while out of the 4 band edges of period 4K, one eigenstate has the form 

ip = dTL~^{x)sn{x)cn^{x) , E = 9 . (48) 

The three other eigenstates of period have the form 

Ip = dn"^(x)cn~^(x)[A + Bsn^ix) + Dsn^{x)] , (49) 

and the corresponding three eigenvalues satisfy the cubic equation 

+ 8(1 + m)r'^ + 80mr + 64m(l + m) = 0, E = -r + 1 + m . (50) 

Further, it is clear that there must be one doubly degenerate state of period 2K whose eigenvalue and 
eigenfunctions are not known analytically. 
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4 GAL Potentials [with at least one parameter a,b,f,g = half-integer] . 

So far, we have discussed GAL potentials when all the four parameters a, b, f, g take integer values. We 
have seen that these are problems with a finite number of band gaps. We now consider the case when at 
least one of the parameters a, 6, f,g is a half-integer. In general, all such problems have an infinite number 
of band gaps and one has only a few QES states. 

4.1 a = half-integer: 

As mentioned in the introduction, we conjecture that the relations © and are simultaneously valid 
when at least one of the four parameters a, b, f,g is a. half-integer. In particular, in case a = k + ^ and 
b, f,g are arbitrary numbers, we assert that k + 1 QES energy values are identical for three GAL potentials, 
that is 

[a = k + ^,b, f,g] 

2(k + b + f + g)+3 2{k + b- f -g)-l 2{k + f -b-g)-l 2{k + g - b - f) - 1 



4 ' 4 ' 4 ' 4 

2{k + b + f-g) + l 2{k + b + g-f) + l 2{k + f + g - b) + 1 2{b + f + g - k) - 1 



■ (51) 



4 ' 4 ' 4 ' 4 

This relation needs some clarification. What is being conjectured here is that there are k + 1 QES mid-band 
energy values of the potential with a = k + 1/2 and 6, /, g being arbitrary numbers, which are the same as 
the band edge energy eigenvalues of the two other potentials given in (|51|) . For these two potentials, we can 
always obtain k + 1 QES band edges, since for both of them the sum of the four parameters characterizing 
the potentials is 2 A; p. 

We have explicitly verified our conjecture in the following cases: (i) Lame potentials with a = 1/2, 3/2 
(and b = f = g = 0); (ii) AL potentials [S] with a = 1/2,3/2, 5 = 1,2,3 (and f = g = 0); (iii) GAL 
potentials with a = 1/2,3/2, and either 6 or / or 17 is arbitrary while the remaining two parameters 
take any integer values. For example, we know that the mid-band state of the a = 1/2 Lame potential 
is at (1 -|- m)/4. Using eq. H51() . we predict that both the potentials [I,— t,— 1,— t] and ^ ^ 



L4' 4> 4> 4J "^^^ L4> 4) 4) 4J 

must have a QES band edge eigenvalue at = (1 -|- m)/4. Using Table 4 of I it is easily checked that 
this is indeed so, with the corresponding band edge eigenfunctions being ip = [dn(x)cn(x)sn(2;)]-^/^ and 
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cn ^/^(j;)sn^/^(x) respectively. We have similarly checked the equivalence in all the above- 
mentioned cases. For example, using the results 16 for the AL potential [3/2, 1,0,0], we find (and verify 
using Table 4 of I) that the GAL potentials [7/4, 3/4, -1/4, -1/4] and [5/4, 5/4, 1/4, -1/4] have QES band 
edges at E = (13 + 5m) /4 it V9 — 9m + m^. The corresponding band edge eigenfunctions are also easily 
written down. Similarly, using the results in I for the mid-band states of the GAL potentials we verify that 
the potentials [t/2 + l,t/2 - 1, (1 - t)/2, -{1 + t)/2] and [(1 + t)/2, -(1 - t)/2, 1 - t/2,t/2] indeed have a 
QES band edge eigenvalue E = t^ + 9m/4, for any non-integer t. 

Thus one is fairly confident that the remarkable relationship (|51j) is indeed valid. We now turn around 
and use H51|) to predict new results for the mid-band states of GAL potentials with half-integer values of 
a in case h,f^g take arbitrary values. As an illustration, we predict that the potential [1/2, 6, /, (7], where 
6, /, g are arbitrary numbers must have a QES state at energy 

S = (6 + l/2)2 + (/ + l/2)2m, (52) 

while the potential [3/2, 6, /, g] with 6, /, g being arbitrary numbers must have two QES states with energies 

= [(6 + 1/2)2 + 1] + [1 + (/ + 1/2) V 
±2^[{h + 1/2) + (/ + l/2)m]2 -{h + f + g + 3/2)(6 + f-g + l/2)m . (53) 

Perhaps some clarification is required as to when the QES state is a band edge and when it is a mid- 
band state. We believe that only those states correspond to band edges for which a + b + f + g [or any 
other combination obtained by changing one or more of these parameters to —a — 1,-6- 1,—/— 1, —g — 1 
respectively] is equal to an even integer (including zero). All other QES states should correspond to mid- 
band states. We thus believe that the QES energy values as given above for the potentials [1/2, /, g\ and 
[3/2, 6, f^g] with arbitrary 6, f^g are in most cases the energies for the mid-band states of these potentials. 
In Sec. 6 we shall obtain the QES eigenstates corresponding to some of these QES eigenvalues by using 
the connection of the GAL potential problem and Heun's equation. 

We can also predict the QES energy values when a = 5/2 and 6, /, g are arbitrary numbers, by computing 
the band edges of either of the two potentials given in 1)51(1 . In this way, we predict that the mid-band 
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energy values for the three QES states of the GAL potential [5/2, b, f, g] are solutions of the cubic equation 

+ 2[1 + 66 + (1 + 6/)m]r2 + 4[2(462 - 1) + 2(4/^ - l)m'^ + (46^ + 4/^ - ig^ + 246/ 
+86 + 8/ - 2^ + 3)m]r + 8m[26 + 1 + (2/ + l)m] [4(6 + ff - {2g + 1)^] = , (54) 

where r = —E + (3/2 — 6)^ + (3/2 — f)'^m. For the special cases (i) 6 = / = 5 = as well as (ii) 
6 = l,/ = g = 0, itis easily checked that eq. ()54() agrees with well known results |TBj, thereby providing 
a powerful check on our calculations. Generalization to higher half-integer values is straightforward (in 
principle) and it is easy to see that energy values for a + 1/2 QES mid-band states can be predicted (at 
least in principle) when 6, /, g are arbitrary numbers and a is a half-integer. 

4.2 a,b = half-integers: 

Let us now discuss the case when both a and 6 are half-integers while / and g are arbitrary numbers using 
our conjecture that eqs. © and ((TJ are both simultaneously valid. In particular, fora = A; + ^, 6 = / + ^ 
while /, g are any numbers, we assert that the three potentials 

1 1 

[a = k + -,b = l + -,f,g] 

k + l + f + g + 2 k + l-f-g k + f -l-g-l k + g-l- f -1 



2 ' 2 ' 2 ' 2 

k + l + f-g + l k + l+g-f + l k + f + g-l l + f + g-k 



(55) 



2 ' 2 ' 2 ' 2 

have the same k + I + 2 QES energy values. This means that the QES mid-band energy values of the 
potential with a = k + \/2,b = 1 + 1/2 are the same as the band edge eigenvalues of the two other 
potentials in ()55() . For these two potentials, we can obtain k -\- I + 2 = a + 6+1 QES states since, for 
both of them, the sum of the four numbers characterizing the potentials is either 2k or 21. This is possible 
because the GAL potential [a,6, /, (/] remains unchanged when any one (or more) of the four parameters 
a, 6, /, g changes to [—a — 1, —6 — 1, — / — 1, —g — 1] respectively. 

As an illustration, consider the potential [3/2, 1/2, /, with /, being arbitrary numbers. It is then 
easily shown that this potential must have 3 QES mid-band states at 

^1 = 4 + (5 + l/2fm , (56) 
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E2,^ = 2 + [(/ + I/2Y + l]m ± J[2 - (2/ + l)m]2 + [{2g + if - (2/ - ifm] 



(57) 



Note that out of the k+l + 2 QES states, the energies for the k+\ states can also be obtained by considering 
the previous case of a = k+1/2 and b, f, g arbitrary and putting b = I + 1/2 at the end of the calculation. 
As an illustration, consider the case [3/2, b, f,g] where b is any arbitrary number. As shown in the last 
section, there are two QES energies given by eq. On putting 6 = 1/2 in eq. H53() we find that the 

two QES energy values of [3/2, 1/2, /, g] are precisely as given by eq. (|57)) . For the special case f = g = 0, 
these eigenvalues agree with well-known results for the AL potential UHI- Similarly, consider the case 
[5/2, b, 0, 0]. As shown in the previous subsection, there are three QES energy values as given by eq. (|54p . 
On putting b = 1/2, it is easily seen that the three eigenvalues are E = (l + 9m/4), (1 + 25m/4), (9 + m/4), 
in agreement with the eigenvalues obtained by us previously |16j . 



4.3 a,b, f = half-integers: 

We shall now discuss the case when three out of the four parameters (say a, b, f) are half-integers while g 
is any number (not a half-integer), and obtain the energies of the QES mid-band states. Our argument is 
again based on the assertion that relations © and Q are simultaneously valid. In particular, we assert 
that fora = /c-|-|,6 = Z-l-|,/ = n-|-| the three potentials 

[a = k + ]^,b = l + ]^J = n+]^,g\ 

2{k + l + n + g) + h 2{k + I - n - g) - 1 2{k + n - I - g) - 1 2{l + n - k - g) - 1 



4 ' 4 ' 4 ' 4 

2{k + l + n-g)+3 2{k + I + g - n) + 1 2{k + n + g - I) + 1 2{l + n + g - k) + 1 



(58) 



4 ' 4 ' 4 ' 4 

have identical k + I + n + 3 QES energy values. What is being asserted here is that the k + I + n + 3 
QES mid-band energies of the potential with half-integer values of a,b, f and arbitrary g are the same as 
the band edge energy eigenvalues of the two other potentials in (|58j) . This happens because for arbitrary 
values of g, for the two potentials 

'2{k + l + n + g) + 5 2{k + I - n - g) - 1 2{k + n - I - g) - 1 2{l + n - k - g) - Y 



4 ' 4 ' 4 ' 4 

2{k + l + n-g)+3 2{k + I + g - n) + 1 2{k + n + g - I) + 1 2{l + n + g - k) + 1 



(59) 
(60) 
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we can always obtain k + l + n + 3 = a + b + f + 3/2 QES band edges, since for both potentials the sum of 
the four numbers characterizing the potentials is either 2k or 21 or 2n. This is possible because the GAL 
potential [a, b, f, g] remains unchanged when any one (or more) of the four parameters a, b, f, g change to 
—a— 1,-6—1, —/ — I, —g — 1 respectively. Thus we conjecture that the potential [k + l/2,l + l/2,n+l/2, g], 
for arbitrary g has k + I + n + 3 QES mid-band states. We now make a number of predictions for the 
energy eigenvalues of the mid-band states for GAL potentials of the form [A; -|- 1/2, Z -|- l/2,n -|- 1/2,^']. 

1. We predict that the potential [1/2, 1/2, 1/2, g] has three QES mid-band states with energy values 

Ei = {l + m) , E2 = l + {g + l/2fm , E3 = {g + l/2f + m . (61) 

2. The potential [3/2, 1/2, 1/2, 5] has four QES mid-band states with energy values 

^1 = (4 + m) , ^2 = (5 + 1/2)^ + 4m , ^3,4 = 2(1 + m)± ^4(1 - + {2g + Ifm . (62) 

3. The potential [3/2,3/2, 1/2, g] has five QES mid-band states with energy values 

Ei = {g + l/2f + 4m , E2,3 = 5 -|- 2m ± ^^4(2 - m)^ -|- {2g + lYm - Am , 

E4^5 = 5 + [{g + l/2f + l]m ± y^[4 - {2g + l)m]2 - {2g - 3)2m + 4m . (63) 

4. The potential [3/2, 3/2, 3/2, 5'] has six QES mid-band states with energy values 

Ei^2 = (5 -h m) ± y^l6(l -m)2 + {2g + lYm , 

EsA = 5 + [ig + 1/2)2 ^ ^ ^16 - {2g + l)^m{l - m) , 

E5,6 = [(g + l/2f + 1] + 5m ± ^16m2 + (2^ + 1)2(1 - m) . (64) 

5. One can readily obtain some QES mid-band energy values when (i) k is an arbitrary integer while I 
and/or n are either or 1 ; (ii) k, I are arbitrary integers with n = or 1 and g being any arbitrary 
number. It would be nice to prove (or disprove) these conjectures and more importantly, try to obtain 
the corresponding energy cigcnfunctions. Wc shall have something to say about this point when we 
discuss the implications of these results in the context of Heun's equation. 
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6. Note that out of the k + I + n + 3 QES states, the energy values for k + I + 2 states can also be 
obtained by considering the previous case [a = k + 1/2, b = I + 1/2, f,g =arbitrary] and putting 
f = n + 1/2 at the end of the calculation. As an illustration, consider the case [3/2, 1/2, f,g] where 
/, g are arbitrary numbers. As shown in the last subsection, there are three QES energy values and 
they are given by eqs. (|56j) and (|57jl . On putting / = 1/2 in these equations, we find that three (out 
of four) QES eigenvalues of [3/2, 1/2, 1/2, g] as given by eq. (|63|1 are correctly obtained. 

4.4 a, b, f,g = half-integers: 

Finally, let us consider the case when all four parameters are half-integers. The potential is of the form 
[k + 1/2, / + 1/2, n + 1/2, p + 1/2], where k, I, n,p are integers and we take k > I > n > p. We shall discuss 
the two cases when the sum A; + / + n + p is an even or an odd integer separately. 

4.4.1 k + I + n + p = even integer: 

As shown in I, for the GAL potential ©, QES energies are obtained when the sum of the four parameters 
is an even integer including zero. In this case, as already shown by Takemura |15j . both the relations © 
and ^ are simultaneously valid. We would like to assert here that in this case, in general there should be 
k + n + l+p + 4: QES states, since the sum of the four numbers characterizing the potentials is either 2k or 
21 or 2n or 2p. This is possible because the GAL potential [a, b, f, g] remains unchanged when any one (or 
more) of the four parameters o, b, f,g change to — o — 1, —b — 1, —/ — I, —g — 1 respectively. Unfortunately, in 
the several specific cases that we have examined, we find that some of the eigenvalues simply get repeated. 
Thus the true number of QES states may be much less than k + l + n + p + 4:. For example, consider 
the case of [5/2, 1, 2, 1/2, 1/2]. While naively we expect six QES states, we only find three QES levels at 
E = (1 + m), (1 + 9m), (9 + m). Similarly, for the potential [9/2, 1/2, 1/2, 1/2] while naively we expect 8 
QES states, we only find four QES levels at = (1 + m), (1 + 25m), (25 + m), 9(1 + m). 
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4.4.2 k + I + n + p = odd integer: 

Let us now discuss perhaps the most intriguing case when all four parameters are half-integers and their 
sum is an odd integer. While it is clear from I that no QES band edges can be obtained when k + l-\-n-\-p is 
an odd integer, it is not obvious whether mid-band states can be obtained in this case. In fact we shall now 
obtain energy values for several QES mid-band states for these potentials. As argued in the introduction, 
while Takemura has shown that relations © or ((T)) are valid when a + b -\- f + g is an even or an odd 
integer respectively (and all are integers), we conjecture that irrespective of whether a + b + f + g is an 
odd or an even integer, both the relations are always valid. Using eqs. © and ((T)), it then follows that 
irrespective of whether a-|-6-|-/ + 5fisan even or an odd integer, they always have a partner potential of 
the form [k + 1/2, 1 + 1/2, n + 1/2, p + 1/2] where k, I, n,p are all integers and their sum is an odd integer. 
We thus conjecture that when a-\-b + f + g is an even integer, then the QES mid-band energy values of 
the potential [{a + b + f + g + l)/2, {a + b - f - g - l)/2, {a + f - b - g - l)/2, {a + g - b - f - l)/2] are 
the same as the band edge eigenvalues of the potential [a, b, f,g]. Similarly, when a + b + f + g is an odd 
integer, then we conjecture that the QES mid-band energy values of the potential [(a + 6 + / — g)/2, (a + 
b + g — f) /2, (a + f + g — b) /2, {b + f + g — a) /2] are the same as the band edge eigenvalues of the potential 
[a,b,f,g]. This is rather remarkable. Since the band edges of the GAL potentials [a,b,f,g] with integer 
a,b, f,g are all known (at least in principle), hence one has a prediction for k + l + n + p + A QES mid-band 
states of the GAL potentials of the type [k + 1/2,1 + 1/2, n + 1/2, p + 1/2] when k + l + n + p is an odd 
integer. To be precise, we predict that the potential [k + 1/2, / + 1/2, n + 1/2, p + 1/2] has k + l + n + p + 4: 
QES mid-band states in case k-\-l + n + p is an odd integer and these eigenvalues are identical to the band 
edges of the potential [{k + 1 + n + p + 3)/2, {k + 1 - n - p - l)/2, {k + n - I - p - l)/2, {k +p - I - n - l)/2]. 
As an illustration, knowing the band edges of the Lame potential 2msn^(x) we predict that the three QES 
mid-band state energy values of the potential [1/2, 1/2, 1/2, —1/2] must be at = m, 1, 1 + m. Similarly, 
we predict that the five mid-band QES energies of the potential [3/2, 1/2, 1/2, 1/2] are at 

E2 = l + m , Es = {1 + 4m) , = {4 + m) , Ei^^ = 2(1 + m) ± 2\Jl -m + ra^ , (65) 

The validity of our conjecture and in a way the consistency of our whole approach can be checked by 
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extrapolating the results obtained in the previous subsection. In particular, in the last subsection we have 
discussed the case when the potential is of the form [k + 1/2,1 + 1/2, n + 1/2, g] where g is any arbitrary 
number and have seen that in that case one obtains energies for k + I + n + 3 QES mid-band states. On 
choosing g = p + 1/2 and choosing k,l,n,p such that their sum is an odd integer, we have verified in many 
cases the validity of our conjectures. 

We thus predict that the potentials [3/2,1/2,1/2,1/2], [5/2,1/2,1/2,-1/2] and [3/2,3/2,1/2,-1/2] 
have exactly the same (five) QES mid-band energies as the band edge energy eigenvalues of the potentials 
[2,0,0,0], [2,1,1,0], [2,1,0,0] respectively. It is worth observing that there are exactly three ways of 
obtaining k + l + n + p= 1 given that k > I > n > p and that while k, I, n are non-negative integers, p is 
> — 1. It may be noted here that in case n is also —1 then the corresponding partner potentials as given 
by eqs. © and ((Tj) are potentials of the type [a, a, 6.6] and similarly if both n = l = —1 then the partner 
potentials are of type [a, a, a, a] which as explained in Sec. 2, are not included in our analysis. Thus, put 
another way, the problem of finding the number of independent potentials with say a band gaps reduces to 
finding four integers k, I, n,p with k > I > n > p (with n > 0,p > —1) such that their sum equals 2a — 3. 

5 GAL Potentials [a,b,f,g = arbitrary numbers]. 

So far, we have discussed the cases when the four parameters a,b,f,g take integer values or at least one 
of them is half-integer. Now we want to extend this discussion to the general case when a, b, f, g take 
arbitrary values. 

As seen in Sec. 4, when either one, two or three of the parameters are half-integer while the remaining 
parameters are arbitrary, then clearly the corresponding partners indeed correspond to the case where 
a, b, f,g are arbitrary numbers. We therefore conjecture that the relations ^ and Q are valid even when 
the four parameters a, 6, /, g take any arbitrary values, the only restriction being that either a + b + f + g or 
any other combination obtained by replacing one or more of these parameters by —a—1, —6—1, —f—1, —g—1 
respectively is a non-negative even integer. Further, in that case the two partner potentials have the same 
band-edge eigenvalues as the QES mid-band energies of the potentials where either one, two or three of 
the parameters are half-integers. 
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A few illustrative examples are in order here. Consider the potential [4/5,2/5,2/5,2/5]. In this case 
while a + b + f + g = 2, no other combination characterizing relations (|38jl to (|46jl gives an even integer. 
Using eq. © we find that this potential has a GAL partner [3/5,3/5,3/5, 1/5]. Using Table 4 of I it is 
easily shown that both these potentials have two (identical) QES band edge energy eigenvalues 

2g 2 

E = —{I + m) ± -Vl - m + . (66) 

ZO 

If instead we consider the potential [17/5, 8/5, 7/5, 6/5], then only a + f — b — g is an integer and the 
corresponding GAL partner potential is [13/5,12/5,11/5,2/5] and both have one (identical) QES energy. 
Of course it can happen that a, b, /, g are such that more than one of the relations H37() to H46|) are satisfied. 
In that case one has more QES band edge eigenvalues. For example, if the potential parameters are such 
that a+b—g — f — 1, a+ f — b—g — 1 as well as a+g — f — b—1 are nonnegative integers, then using eqs. to 
H4U|) it is easily seen that the number of QES energy eigenvalues is equal to {3a—b—f—g) /2. One illustration 
of this is the potential [11/5, 1/5, 1/5, 1/5] which has three QES energies (1+m), 144/25+m, l + (144/25)m. 
We might add here that the discussion above is valid even if the numbers a, b, /, g are irrational numbers 
but such that a + b + f + g or any other combination obtained by replacing one or more of a,b, f,g to 
—a — 1, —b — 1, — / — 1, —g — 1 respectively is an even integer (including zero). 

We would like to restate here that when all four parameters a, b, f, g are integers, one has a finite 
number of band gaps. In all other cases one expects to have an infinite number of bands and band gaps 
out of which only a few are QES states. 

6 Implications for Heun's Equation. 

Heun's equation, a second order linear differential equation with four regular singular points has been 
extensively discussed in the mathematics literature |2U1 1211 1^ . The intimate connection between Heun's 
equation and GAL potentials is well known j^I. In recent years, this equation has also proved very useful in 
the context of a number of physical problems, like quasi-exactly solvable systems |23| . sphaleron stability 
j24j . Calogero-Sutherland models H^j, higher dimensional correlated systems |26| . Kerr-de Sitter black 
holes ^7\, and finite lattice Bethe ansatz systems j28j . 
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The canonical form of Heun's equation is given by [20] 



G{x) = , (67) 



dx'^ X x — 1 X — c dx x{x — l){x — c)_ 
where a, (3, 7, 6, e, q, c are parameters, except that c / 0, 1 and the first five parameters are related by 

7 + 5 + e = a + /3 + l . 

The four regular singular points of eq. (|67|) are located at x = 0, 1, c and the point at infinity. 
If we make the transformation x = sn^ (y,m), then Heun's equation takes the form j2Uj 

F"{y) + [(1 - + (1 _ ,^)Sn(y,m)dn(y,m) ^ _ ^^cn{y,m)dniy,m) 

dn(y,m) cn(y,m) sn(y,m) 

— [Amq — 4a/5msn^(y, m)]F(2/) = , (69) 

where c = 1/m, G{x) = F{y). The periodic solutions of eq. (|69() correspond to the polynomial solutions 
of eq. 1)6 7|) while the quasi-periodic solutions correspond to non-polynomial solutions of (|67() . 

The interesting point is that after a transformation, the Schrodinger equation Q for the GAL potential 
((21) is in fact Heun's eq. In particular, let us start from the Schrodinger equation @ for the GAL 

potential ((3). On substituting 

V'(y) = dn-^(y)cn-^(y)sn-9(y)0(y) , (70) 

one can show that (j){y) satisfies the differential equation 

\ , or , sn(y,m)cn(y,m) sn(y, m)dn(y, m) cn(y,m)dn(y,m) , 

(f) [y) + 2[mb — ^ + f ^ g ^ J0 (y) 

dn(y,mj cn[y,mj sn(y,mj 



[i?-Qmsn^(y,m)]0(y) =0 , (71) 



where 



R=-E + m{g + bf + {f + gf, Q = (b + f + g){b + f + g - I) - a{a + 1) . (72) 

Thus once we obtain solutions of the Schrodinger equation for the GAL potential ^ , then we can imme- 
diately write the solutions for the periodic form of Heun's eq. (|69jl and the solutions of the original Heun's 
eq. (|67|) with the identification 

l = --g,6=--f, e--b, 

a + p=^-{b + f + g) , 4a(3 = Q , Amq = R ,F{y) = <P{y) . (73) 
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We now make a crucial observation. From eq. H72j) . it follows that if under any transformation, the 
parameters bi, fi,gi change to b2, /2;52 and the energy E remains invariant, then the corresponding values 
of R are related by 

i?i - m(6i + gif - (A + gif = R2 - m{b2 + 92? - U2 + 92? ■ (74) 

Making use of eq. H74|) and the connection between GAL potentials as given by eqs. © and ((T)), we can 
obtain interesting relations for Heun's equation. Using the fact that the two GAL potentials given by eq. 
(jEJ have the same band edge energy eigenvalues and the eigenfunctions for both the partners as are given 
by eqs. 1)3 7|) to H41|) . one can obtain the connection between the two corresponding solutions of Heun's 
equation. For example, consider the solution ()37|) and the corresponding solution of the GAL partner 
obtained by the above substitution. Using eq. H74|) it then follows that corresponding to a given periodic 
(i.e. polynomial) solution of Heun's equation with parameter set (a, /5, 5, e, 7, q) there always exists another 
periodic solution with the same q provided the other parameters change as follows: 

7^a, a^7, P ^ P , e^l + /3-5,(5^1 + /3-e. (75) 

Similarly, on considering the other three periodic solutions as given by eqs. ()38|) to (|4U() and the corre- 
sponding solutions of the GAL partner potential with the same energy, we find that corresponding to a 
given periodic solution of Heun's equation, there exist the following three (periodic) solutions with the 
change of parameters given by (note that R = 4mq and c = 1/m) 

^q-l3{5-a) (76) 

q-f3{e- a)c , (77) 

q + (3[{a-6) + ia-e)c] . (78) 

Thus given a periodic solution of Heun's equation, one immediately has 4 other periodic solutions as 
given by eqs. (|75() to H78|) . We have checked that if instead we consider the two partner GAL potentials 
given by eq. ((T)) and consider the corresponding eigenfunctions given by eqs. 1)42(1 to ()46() (and those of 
the corresponding GAL partner potentials with the same energy), then we again obtain the same relations 



7^1 + /? — e, a— >5, f3 ^ (3 , e^l + /3 — 7,5— >a, 



e ^ a , a^e, (3 ^ fJ , 7^1 + /3 — (5,(5^1 + /3 — 7, q 



7^1 + /3-7, a ^ 1 + P-a , P ^ P , e^l + P-e,6^1 + P-S,i 
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[eqs. (|75)) to ((75|) ]. As an additional check, we have looked at the partner GAL potentials as given by eqs. 
(|51|) . (|59jl . (|6Ujl . (|55|) and © and in all these cases we get back the relations (|75|) to (|78|). which to the best 
of our knowledge, are new results. Several comments are in order: 

1. In Sec. 4 we have shown that the QES mid-band energy values of the GAL potential [a = k + 
1/2,6, 7,(7] are the same as the QES energies of the two GAL potentials given in What does 

this imply in the context of Heun's equation? It is easily shown that as a consequence of the discussion 
in Sec. 4, given a periodic solution of Heun's equation with the set of parameters a, /3, 7, 6, e, q, one 
has a corresponding quasi-periodic solution with changed parameters: 



An additional check on this relation is obtained by using the connection between the mid-band states 
of the GAL potentials [a = k + 1/2, b = I + 1/2, f,g] and [a = k + 1/2, b = I + 1/2, f = n + l/2,g] 
and the QES energies of the potentials and (I^Uj) and we again obtain the same connection 

between the periodic and quasi-periodic solutions of Heun's equation. 

2. Using the results in Sec. 4 regarding the case when three of the four parameters a,b,f,g are half- 
integers, we also obtain two more relations connecting the periodic and quasi-periodic solutions of 
Heun's equation. In particular, given a periodic solution of Heun's equation with the set of parameters 
a, (3,'y,6,e,q, it implies the following two quasi-periodic solutions of Heun's equation: 

7^a, a ^ 1 + a - e , (5^5, e^l + 5- l3,5^1 + l3-€,q^q + a{(3 - 5) , (80) 

7^a, a^l + a-6 , f3 ^ e , e^l + P-6,6^1 + e-P,q^q + 4a(/3 - e)c . (81) 

3. Needless to say that if instead, a quasi-periodic solution of Heun's equation is given, then by inverting 
eqs. ()79() to (|81() . we immediately obtain three periodic solutions of Heun's equation. 

So far we have discussed how the connections between different GAL potentials can help in finding new 
solutions of Heun's equation. It may happen that in some cases it may be simpler to solve the algebraic 



7^2 — Q, a^l-l-7 — a, f3 —>■ 1 + P — a , e 



l + P-5 ,5^1 + P-e 



q^q+{a-l)[{l+/3-6) + {l + /3- e)c] . 



(79) 
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Heun's equation (|67() rather than its periodic variant. We now show that this is indeed so in the case of 
several quasi-periodic mid-band eigenfunctions. Consider for example the GAL potentials when either b 
or / or (7 is 1/2 while the other three parameters are arbitrary. Using arguments of section IV, we can 
easily obtain the eigenvalues for mid-band states for these potentials. We shall now show that using these 
eigenvalues we can easily solve the algebraic form of Heun's eq. (|67|) and hence using the connection as 
explained above, obtain the eigenfunctions for the mid-band states of these GAL potentials. 

Consider the GAL potential [a,l/2, f, g] where a,f,g are arbitrary numbers. Note that using the 
relations Q Table 4 of I, it is easily shown that the QES mid-band eigenvalue of the GAL 
potential [a, 1/2, f,g] is at ^ = {a + l/2f + m{g + 1/2)^ . On using the connection formulas H72|) and 1)73(1 
it is easily shown that the corresponding parameters for Heun's eq. ()67|1 are 

7 = 1/2-5, 5 = 1/2-/, e = 0,a = (a - 5 + l/2)/2 , 
/? = -(a + / + g + l/2)/2 , g = (a + / + 5 + l/2)(/ + 5 - a - l/2)(c/4) . (82) 

Remarkably, for these parameters, it is straightforward to obtain the solution of the algebraic Heun's eq. 
(|S7|) and show that 

G{x) = F[{a -f-g + l/2)/2, -{a + f + g + l/2)/2, 1/2 -g;x], (83) 

where F(a, b, c; x) is the hypergeometric function. The corresponding mid-band state eigenfunction for the 
GAL potential [a,l/2, f, g] is then immediately written down. We have verified that this is indeed the 
correct eigenfunction in the following cases (i) a = / = g = (ii) a integral, f = g = (iii) a arbitrary 
while f,g are integral. 

Similarly, for the GAL potential [a, b, 1/2, g], with arbitrary a, b,g, the QES mid-band energy eigenvalue 
is E = {g-\-l/2)'^ -\- (a-|-l/2)^m and proceeding as above, it is easily shown that the solution of the algebraic 
Heun's eq. (|H7)) is given by 

G{x) = F[{a-b- g + l/2)/2, -{a + b + g + l/2)/2, 1/2- g;mx] . (84) 

And finally, for the GAL potential [a,b, f,l/2], with arbitrary a,b,f, the QES energy is E = (/ + 
1/2)^ + {b + 1/2)^771 and proceeding as above, it is easily shown that the solution of the algebraic Heun's 
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eq. (|H7|) is given by 

G{x) = F[{a -b-f + l/2)/2, -{a + b + f + l/2)/2, 1/2 - 6; (1 - mx)/{l - m)] . (85) 

On using the solutions (|83|) to (jHSJ) it follows that for the potential [a, 1/2, 1/2, 1/2] one knows three QES 
mid-band energy eigenstates. In the special case when a = 2k+3/2 these eigenstates are the mid-band QES 
eigenstates for the potential [2k + 3/2, 1/2, 1/2, 1/2] where the sum of the four parameters characterizing 
the potential is an odd integer. 

7 Summary and Open Questions. 

In this paper, we have addressed many issues regarding GAL potentials with a number of choices for the 
parameters a, b, /, g. The most interesting case is when all the four parameters are integers. This is a 
potential with a finite number of band gaps. We have been able to count the number of independent 
GAL potentials with a given number of band gaps and completely specify the nature of the band edge 
eigenf unctions. We have introduced the new concept of self-dual potentials which are not self-isospectral. 
We are also able to specify how many of the independent potentials with a given number of band gaps have 
supersymmetric partner potentials and how many have non-supersymmetric partner potentials. Finally, 
using the results for the GAL potentials, we have shown that given any one periodic solution of Heun's 
equation, one can obtain four more periodic solutions. 

We have also discussed several issues related with GAL potentials when one or more of the pa- 
rameters take half-integer values. In particular, while nothing is known so far about GAL potentials 
when three of the parameters take half-integer values, we have been able to obtain the QES energy val- 
ues for several of these potentials. Further, using these eigenvalues and the algebraic form of Heun's 
equation, we have also been able to obtain the corresponding eigenfunctions for potentials of the form 
[a, 1/2, /, g], [a, b, 1/2, g], [a, b, /, 1/2] where a, b, /, g are arbitrary numbers. The key point to make while 
addressing these questions is that the relations © and ((T)) are not only valid when the four parameters 
a, b, f,g are integers but also when one or more of these parameters take half-integer values. This in turn 
immediately implies that these relations are also valid when the four parameters a, b, /, g take arbitrary 
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values so long as either their sum a+b+f+g [or one or more of the combinations obtained by changing one 
or more of the parameters to —a — 1,—b — 1, — / — 1, —g — 1 respectively] is a nonnegative even integer. We 
have also conjectured that both relations © and Q are simultaneously valid when a, b, /, g are integers 
and that the energy eigenvalues for the band edges of these potentials are the same as mid-band QES 
energy values of GAL potentials in which all four parameters are half-integers and their sum is an odd 
integer. Finally, using these results we have also shown that given a periodic solution of Heun's equation, 
one can immediately obtain three quasi-periodic solutions of the same equation. 

This work raises several issues which we have not been able to address satisfactorily: 

1. Can one explicitly write down all seven KdV equations of seventh order? 

2. What are the QES eigenfunctions for GAL potentials when one or more of the parameters is half- 
integral (> 3/2) while the remaining parameters are arbitrary? 

3. The problem when two of the four parameters are half-integers needs further study. In particular, it 
is still not clear how many QES energy eigenvalues can be obtained, in general, in that case. 

4. When the sum of all the four parameters is an even integer, it is clear that the QES states correspond 
to band edges. However, a complete understanding is still lacking regarding the number of QES states 
for various values of a,b,f,g. Further, when the sum of the four parameters is an odd integer, the 
form of the QES eigenfunctions is not clear when the half- integer parameters are > 1/2. 
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